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Abstract 

We describe the structure of bifurcations in the unbounded classical 
Diamagnetic Kepler problem. We conjecture that this system does not 
have any stable orbits and that the non-wandering set is described by a 
complete trinary symbolic dynamics for scaled energies larger then = 
0.328782 .... 

The problem of a charged particle moving in a combined attracting electri¬ 
cal Coloumb field and a strong homogeneous magnetic field is a simple system 
showing complicated chaotic behavior. The corresponding quantum mechanical 
problem is a Hydrogen atom in a strong magnetic field. This system is called 
the Diamagnetic Kepler problem (DKP) and in the last decade this problem 
has been much investigated both classicly and quantum mechanically in theo¬ 
retical and experimental studies. Review articles by Friedrich and Wintgen[l] 
and Hasegawa et.al.[ 2 ] yield good introductions to this problem. One interest¬ 
ing result of this research is an astonishing relation between the classical and 
the quantum mechanical systems, where the quantum mechanical energy lev¬ 
els analyzed in a proper way reflect the periodic orbits existing in the classical 
system. Most of these results have been on the bounded system with negative 
energy such that the electron cannot escape. We report here some numerical 
results concerning the classical orbits in the DKP with a positive energy, where 
quantum mechanically the atom can ionize. We also give some remarks on the 
implications for the quantum mechanical systems from the results obtained here. 

The Diamagnetic Kepler problem with zero angular momentum is given by 
the Hamiltonian 

H = I2me - /r+ ^meLo‘^{x^+y‘^) ( 1 ) 

Using suitably scaled variables and semi-parabolic coordinates this problem can 
be turned into the study of the Hamiltonian [1, 2] 

h=]^pl + + p^)= 2 ( 2 ) 

at the fixed “pseudo energy” 2 and with a scaled energy e = depending 

on both the energy, E, and the strength of the magnetic field, 7 . We consider the 
scaled energy e to be a parameter when we investigate the classical trajectories. 

The Hamiltonian (2) has a rather complicated bifurcation structure for e < 0 
and there is not any known symbolic dynamics description of the orbits. For 
e > 0 it has been proposed by Eckhardt and Wintgen[3] that the orbits may 
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be described by the same symbols as the 4-disk system [4] and that this 4- 
disk symbolic dynamics for periodic orbits (cycles) can be obtained in DKP 
by finding self conjugate points. The 4-disk system consists of a point particle 
moving freely on a 2-dim. plane and bouncing elastically off the boundaries of 
4 disks in the plane. The 4 disks have radius 1 places with the disc centers at 
the four points (±r/2,±r/2) int the (i/,/u)-plane. The center-center distance r 
is the parameter for the 4-disk problem. 

For the Hamiltonian system with the smooth potential V{x,y) = it 

has been argued that the bifurcations are very similar in this smooth potential 
and in the 4 disk billiard [7, 8]. This seems to be true also for the DKP with 
positive e and this yields a method to determine stable orbits and to find the 
parameter value where the repellor is a complete Cantor set. We can then 
show that there are no elliptic orbits above a critical energy Cc which has to be 
determined numerically. 

In a dispersing billiard system the singular orbits which define the pruning 
front [9, 10, 11] are the orbits tangent to the boundary. The pruning front 
determines which orbits that are admissible in the system and all bifurcations 
take place on this pruning front. In the 4-disk billiard we obtain a symbolic 
dynamics description by enumerating the disks clockwise and the symbolic dy¬ 
namics of a trajectory is given by a string ... s_iSoSiS 2 • • • where St is the label 
of the disk the particle bounces off at the integer time t. Other alphabets 
can be defined from this basic alphabet [4, 11]. The first orbits which become 
forbidden in the 4-disk system as the distance between the disk centers de¬ 
creases to a critical distance Tc = 2.2046..., are the heteroclinic orbits 14 21 
and 14121 (a line above a string denotes an infinite repetition of this string 
14 21 = ••• 141414212121 •• •) and orbits symmetric to these [10, 11]. These 
orbits are the corners of the 2-dimensional Cantor set of non-wandering[5] (not 
escaping) orbits in the phase space and are the first to become non-admissible 
in a symmetric system. We show below that these heteroclinic orbits also are 
the corners of the non-wandering Cantor set for the complete chaotic DKP. 

The idea that singular orbits determine the set of admissible orbits can be 
carried over from the hard billiard systems to the soft potential systems. One 
then has to determine the orbits in the smooth system which correspond to the 
tangential orbits in the billiard. A second way to view the analogy between hard 
and soft systems is to draw the manifold structures in a Poincare map. This 
yields a simpler and maybe more familiar picture, but the symbolic dynamics 
and physical interpretation are more difficult to get from this point of view. We 
will here discuss both approaches to the problem. 

There does not exist a steep wall in the smooth system as it does in the 
billiard system and consequently no direct way to find a “tangent” orbit, but 
the singular orbits in the billiards have another property to which we can find 
an analog to in smooth systems and we will apply this to determine the singular 
orbits. Assume that there exists a singular orbit tangent to the boundary in 
the 4 disk system. Then choose a second orbit close to the singular orbit at 
the tangent point but a finite distance S away from the boundary and such 
that this second orbit does bounce off different disks than the first orbit after 
some time. This second orbit will then have a different symbolic description 
from the singular orbit. Then there will also exist a third orbit bouncing off 
the boundary near the tangent point of the first orbit with some finite angle ^ 
with a symbolic description identical to the second orbit both in the future and 
in the past. We can define the singular orbit as an orbit with extremum future 
and past trajectories since the two different neighbor trajectories have identical 
symbolic dynamics in past and future. In a well ordered symbolic alphabet [11] 
this yields a extremum symbolic value and in the configuration space the two 
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neighbor orbits moves on the same side of the singular trajectory. This new 
definition of singular orbits as extremum orbits can be applied also in the case 
of a smooth potential. 

The set of all singular orbits is a one-parameter family of orbits in the 4-disk 
system and this seems to be true also for the smooth Hamiltonian (2) in the 
DKP. 

In the following we for simplicity mostly discuss the singular orbit normal 
to the diagonal = v. It is clear from the symmetry of the problem that there 
will exist such a singular orbit and as for the 4-disk system this singular orbit 
becomes the heteroclinic orbit 14 21 at the critical parameter value ec. Then 
there is no pruning for e > Cc- The singular orbit on the r'o = fJo line can be 
found numerically by scanning a one dimensional starting point r'o = /J-o keeping 
Pvo = ~Pmo while for other non-symmetric singular orbits one has to scan both 
the starting points and the direction of the velocity. 

Fig. 1 shows a segment of the singular orbit for the scaled energy e = 0 from 
the starting point = Po = 0.816 together with two neighboring orbits starting 
parallel with this at vq = Po = 0.716 (dashed curve) and at vq = po = 0.916 
(dotted curve). Both the two neighbor trajectories are on the same side of 
the singular trajectory after some transient time both in the future and in the 
past. There are no trajectories on the other (the forbidden) side of the singular 
trajectory both in the past and in the future. There may be orbits on this 
forbidden side in the past or in the future, but not both in the past and in the 
future. 

The orbit starting at vq = po = 0.716 corresponds to the orbit in the billiard 
system which passes the boundary of disk no. 1 at some distance S, while the 
orbit starting at Uq = po = 0.916 corresponds to the billiard orbit bouncing 
off the boundary of disk no. 1 with some angle (ft. By drawing the orbits in 
configuration space we can determine the singular orbit and we find that the 
singular orbit for e = 0 starts at vq = Po = 0.8166179_ 

The singular orbit in Fig. 1 shows that there cannot exist any orbit described 
by a symbol string containing a substring _414212_ or _4141212_ for e = 0. This 
orbit would be on the forbidden side of the singular trajectory. One forbidden 
periodic orbit is then 41412121. Since the symmetric singular orbit for e = 0 
numerically appears to be non-periodic there is an infinite number of forbidden 
symmetric sub-strings of increasing length determined by this singular orbit. 

To determine all forbidden sub-strings one also has to determine the non- 
symmetric singular orbits. One example of a non-symmetric singular orbit is 
drawn in Fig. 2. We have here numerically determined a trajectory with ex¬ 
tremum future and past compared to orbits starting parallel to it. This orbit 
starts at vq = 0.56166, po = 1.56166 and with Pv^IPixq = — tan(0.605). From 
Fig. 2 we can read off that for example the orbits containing a substing _1421212_ 
or _14121212_ are forbidden. This is a slightly different rule from that we got 
from Fig. 1. All the rules obtained from all singular orbits yield all rules which 
determine the forbidden orbits in the system. 

As e increases the singular orbits will change and for some values of e a 
singular orbit may return close to the starting point with a similar direction. 
When the singular orbit returns it implies that we get a stable periodic orbit. 
Since orbits near the singular orbit do not diverge but approach each other 
around the starting point of the singular orbit this gives rise to stability for 
periodic orbits. One example of this is given in Fig. 3 for the scaled energy 
e = 0.1113 where the singular orbit is found at po = 0.8617. The stable cycle 
existing here is 41412121. This cycle we found above to be inadmissible at e = 0 
and it will be admissible for all e > 0.1113. The cycle is in the center of a stable 
island and it is not directly a forbidden and legal side of the singular orbit here. 
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but the admissibility is determined by the rotation number of the elliptic orbits. 
The nearby orbits drawn in Fig. 3 start outside the stable island. 

Inserted in Fig. 3 is a Poincare map showing when the trajectory 

crosses the v = p, line for 5 different starting points close to the singular tra¬ 
jectory (0.855 < V < 0.869,1.4354 < < 1.4375). We find here the usual 

Hamiltonian structure of KAM-curves and resonances. We conjecture that all 
the periodic and quasiperiodic orbits within the stable island can be given a 
unique symbolic coding by adiabatically following the trajectories into the com¬ 
pletely chaotic region e > ec, and that they belong to one bifurcation family ot 
the 4-disk system, se ref. [8]. The symbolic dynamics of these orbits are given 
as a string • • • 414s_i212so414si212 • • • where St is the symbol 1 or no symbol. 

All the cycles symmetric across the p = v line are stable when they are on 
the symmetric singular orbit and hyperbolic otherwise. 

We will now try to find the parameter value ec where the singular orbit 
escapes directly to z/ oo in the future direction and to /z oo in the past 
direction without ever returning. It is however not trivial to decide when a 
trajectory leaves to infinity in this system because for an arbitrary large value 
of v (or of /i) there exist orbits coming from and returning to the v = 0 line 
(/i = 0 line). One example of this is an orbit crossing the zz-axis with p^ = Q for 
a large value of v. This trajectory will cross the zz = 0 line both in the future 
and in the past and may be a not escaping trajectory. 

Numerically it turns out that it is possible to distinguish escaping and not 
escaping trajectories by plotting the Poincare map (/U, p^) for each crossing of the 
p-axis with > 0 (correspondingly for v). Escaping trajectories yield points 
approximately along or above a straight line p^ = afx + b with a > 0. Returning 
trajectories will yield points along a oval shaped curve in the Poincare map 
(p,p^). We can then distinguish between escaping and not escaping trajectories 
by following the trajectory a reasonably short time. 

Fig. 4 shows the Poincare map for three orbits starting at p = z/ = 0.75 
(O), 0.9 (-I-) and 1.05 (□), p^^ = —p^o) the scaled energy e = 0.335. The 
trajectory from 0.9 escapes directly while the other two return to the p = 0- 
line (but may escape at some later time). The escaping trajectory is drawn in 
configuration space (zz, p) in Fig. 5. The orbit simply comes in from p = oo and 
continues out to z/ = oo. When the symmetric singular orbit escapes directly 
there is no pruning left in the system and we will have a complete Cantor set 
repellor. Other non-symmetric singular orbits will for some e < ec escape to 
infinity either in the future or in the past. There will not be any stable cycles 
since stable cycles only appear near a non-escaping singular orbit. 

There is one orbit missing from the Cantor set and this is the cycle 12 which 
is only approach ed in the v ^ oo limit. Orbits approaching this cycle, for 
example 34(12)"; n oo, are then not bounded by a maximum value of v, but 
they are unstable. 

Numerically we determine the critical parameter value to be ec = 0.328782 .... 

A different way to discuss this problem is to draw the manifold structures in 
a Poincare plane. In this way we can visualize how the stable and the unstable 
manifolds cross each other with a finite angle for e > Cc. When all stable and 
unstable manifolds cross with a finite angle (no tangencies), then the system 
is hyperbolic. In Fig. 6 (a) we have drawn all starting points in a Poincare 
plane not escaped after some fixed integration time T for the scaled energy 
e = 0.35. We have chosen the Poincare plane where the trajectory crosses the 
zz = p line from z/ < p to z/ > p with the coordinates (p, ^), where (p is the 
angle of the trajectory with the zz-axis; Pfi/pi, = tanz/i. This Poincare plane is 
finite with —37r/4 < 0 < 7r/4 and —Pmax < p < Pmax where Pmax(e) is found 
analytically from (2). To determine if a trajectory has escaped we use the idea 
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from Fig. 4 that an orbit escapes if the value of (p^) increases from one to 
the next crossing of the p-axis (z^-axis). The stable manifolds, IT®, are included 
in the set of not escaping trajectories, and for a sufficiently large T this yields 
a good picture of TT®. All points on the border curve, (f) = —Srr/d, 0 = 7r/4, or 
p = ±/ijf,ax in the Poincare map, correspond to the hyperbolic period two cycle 
13. The structure close to the border can be mapped into the structure in the 
middle of the Poincare plane (close to the period two cycle 24) by e.g. v —v. 
Below we also draw the manifolds in a Poincare map that does not have this 
strange border but which on the other hand is infinite. 

The unstable manifolds TT“ are symmetric to IT® with respect to the line 
(p = —7r/4. Both and IT® are drawn in Fig. 6 (b) and the crossing points 
between the and IT® curves yield the non-wandering set. We find here 
that all the crossings between the curves of and IT® have a finite angle and 
consequently: all orbits are hyperbolic. In Fig. 6 it is simple to recognize a three- 
interval Cantor set structure of the manifolds. It is slightly more complicated to 
recognice this at the edges because the border curve corresponds to one unstable 
cycle, and we have here only one half of the Cantor set. This only depends on 
our choice of Poincare map and does not give any real problems. 

These figures give a similar structure of manifolds as the manifold structure 
one can find for a 4-disk system. The stable and the unstable manifolds for the 
4-disk system in a corresponding Poincare map are drawn in Fig. 7. In Fig. 7 (a) 
both the stable and the unstable manifolds are drawn for the disk-disk distance 
r = 2.35, —37r/4 < (f> < 7r/4, l/\/2 — r/2 < p < r/2 — l/v^. This yields nice 
trinary Cantor sets as crossing points between IT“ and IT® in the same way 
as we found it in Fig. 6 (b). The main structure is exactly the same in the 
billiard system. Fig. 7 (a), and in the smooth system. Fig. 6 (b). The important 
difference is that in DKP all folds smoothly connect for large values of \p\ while 
for the 4-disk system the manifolds are discontinuous at the line p = Pmn^- 
This discontinuous jump from ^ to 7r/2 — ^ for the folds of the 4-disk system is 
loosely speaking for DKP glued together and the folds smoothly connect at the 
line 6 = —7r/4. The two drawings yield the same number of crossings between 
IT“ and IT®. 

We can also compare the pruned systems for the DKP (e < tc) and the 
4-disk system (r < Vc) in this Poincare plane. In Fig. 7 (b) the stable manifold 
IT® for the 4-disk system is drawn for the disk distance r = 2.06. We have 
here drawn a quarter of the full phase space —7r/4 < (p < 7r/4, 0 < p < 
r 12 — l/v^- We find that the largest white area which corresponds to orbits 
escaping without ever returning, does not any more reach the p) = —7r/4 line 
as it does in Fig. 7 (a). Some crossings between the stable and the unstable 
manifolds existing for r = 2.35 are here lost, and then also an infinite number 
of periodic orbits are inadmissible. A sketch of the manifold structure are given 
in Fig. 8 (a). The points where a cusp of the stable manifold touch a fold of the 
unstable manifold are marked with circles. These singular points are exactly the 
images of the tangent orbits. Moving the disks slightly closer to each other makes 
the homoclinic orbit at this touching point inadmissible. This is the analog in a 
billiard to the point called a homoclinic tangency for smooth dynamical systems. 
For simplicity we also denote this point a homoclinic tangency. One important 
remark here is that the homoclinic tangencies indicated in Fig. 8 (a) are the 
primary tangencies and all other tangencies in other part of the phase space are 
images of these. The different primary tangencies are associated with different 
tangent orbits in the billiard. There is a one to one correspondence between a 
primary homoclinic tangency and singular orbits in the billiard. The homoclinic 
tangencies can be used to define the pruning front. 

In Fig. 9 manifolds are drawn for the Diamagnetic Kepler problem with 
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e = 0.2 < ec- In Fig. 9 (a) both the stable and the unstable manifolds are 
drawn while in Fig. 9 (b) only the stable manifolds are drawn in one quarter 
of the phase space 0 < /i < /imax, —7r/4 < (/> < 7r/4. Comparing Fig. 9 (b) 
with Fig. 7 (b) we find a very similar structure. The main difference is that the 
mainifolds in Fig. 9 (b) are smooth while in Fig. 7 (b) the bending points are 
sharp. For the DKP we have homoclinic tangencies where the unstable and the 
stable manifolds are tangents and for an infinitesimal smaller parameter e these 
homoclinic orbits are inadmissible. The structure of the homoclinic tangencies 
is sketched in Fig. 8 (b). As for the 4-disk system we have a one parameter 
family of primary homoclinic tangencies each associated with a singular orbit. 
The situation in the soft potential is slightly more complicated because a stable 
orbit creates a stable island in the neighborhood of a homoclinic tangency. 
One example of such an island is the orbit in Fig. 3. One way to numerically 
determine the pruning front is to draw the unstable and the stable manifolds 
and determine the set of primary homoclinic tangencies. This can also be done 
for the closed system (e = 0) by using the manifolds of a hyperbolic periodic 
orbit. 

To further illustrate the hyperbolicity for this system, we draw the manifolds 
in a different Poincare map where the trajectory crosses the /i-axix with v 
changing from negative to positive. The manifolds in Fig. 9 (b) for the DKP 
at e = 0.20 are drawn in Fig. 10 (a). Here there is a number of homoclinic 
tangencies giving rise to pruning and stable orbits. The manifolds in Fig. 6 (b) 
for the hyperbolic map at e = 0.35 are drawn in Fig. 10 (b). Here we find that 
the manifolds always cross with a finite angle and there are no tangencies. In 
Fig. 10 (c) the same manifolds are drawn for a large scaled energy, e = 2.0. 
Here we find that the same Cantor set of crossings between the stable and the 
unstable manifolds but the Cantor set is thinner. It is not created any new 
non-escaping orbits when e increases because the manifolds only moves closer 
to the (p = 0 and (p = tt lines. It is reasonable to conjecture that no new orbits 
are created as e increases further and the non-wandering set of the DKP has 
the same topological structure for all €> 

An advantage of this Poincare map is that we avoid the complicated border 
orbit we had in Figs. 6 and 9, and it is easy to see that no new orbits are 
created as e increases. The disadvantage is that the plane is unbounded and 
as II increases we get more and more copies of the foliations. Each new copy 
only has one more bounce than the presiding one and does not yield any new 
information. 

In Fig. 11 equipotential lines are drawn at tc- In addition to the hyperbola 
like curves also existing at e = 0, there are saddle points and a “hill” at the 
origin. These latter structures turn out not to give rise to any new singular 
orbits and stable cycles. The potential at the origin is 0 and trajectories cannot 
turn there. The hill structure only adds instability to the Cantor set repellor. 

In quantum mechanical calculations [12, 13, 14] one has found resonances 
and some with a very narrow width. In quantum mechanical calculations these 
appear because of cancellation of some matrix elements. Semiclassically these 
can be associated with stable orbits like the one in Fig. 3. When we have a pure 
hyperbolic classical system (e > Cc) we can predict that a quantum mechanical 
cross section (for e constant) does not have any narrow width resonances. All the 
other resonances can be calculated by periodic orbit expansion of the unstable 
periodic orbits and work on this continues. 

The DKP with negative scaled energy has a more complicated structure 
and more mechanisms giving stable orbits. For small negative energies the 
pruning front orbits discussed above still exist as a part of the complete pruning 
mechanism. 
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From the results presented above it seems like the Hamiltonian (2) with 
e > ec is a completely chaotic scattering smooth system, without any singu¬ 
larities and where the set of all periodic orbits (and the non-wandering set) is 
unbounded in the configuration space. To the authors knowledge there are no 
earlier examples of this in the literature. 

The author wish to express his gratitude to Dieter Wintgen for the few 
interesting discussions time allowed. The author thanks friends in Copenhagen 
and Freiburg for good discussions. The author is grateful to the the Alexander 
von Humboldt foundation for financial support. 
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Figure 3: A singular orbit and two neighbor orbits for e = 0.1113 where the singular 
orbit is a stable periodic orbit. Inserted is a Poincare map (/x, </<) for v = showing 
KAM-curves. 



Figure 4: The Poincare map {pi,Pv) for v = Q of three trajectories from i/q = 
Mo,Piyo = —Pno showing that there is one trajectory escaping directly yielding a 
completely chaotic system. 
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Figure 6: The manifolds of DKP in the Poincare map {fi, (ft) for v = for e = 0.35. 
(a) TF*, (b) IF* and TF“. 
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Figure 9: The manifolds of DKP for e = 0.2. (a) IF® and fF“ (b) IF®. 
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Figure 10: The manifolds of DKP TF^Hnd in the Poincare map for 

jy = 0 for (a) e = 0.2, (b) e = 0.35 and (c) e = 2.0 
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